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We discuss the extension of the superflavor symmetry of doubly heavy baryons to states which 
contain an excited heavy diquark, and we examine some of the consequences of this symmetry for 
the spectra of doubly heavy baryons and heavy mesons. We explore the ramifications of a proposed 
symmetry that relates heavy diquarks to doubly heavy mesons. We present a method for determining 
how the excitation energy of a system containing two heavy quarks will scale as one changes the 
strength of the interactions and the reduced mass of the system. We use this to derive consequences 
of the heavy diquark-doubly heavy meson symmetry. We compare these consequences to the results 
of a quark model as well as the experimental data for doubly and singly heavy mesons. We also 
discuss the possibility of treating the strange quark as a heavy quark and apply the ideas developed 
here to strange hadrons. 

I. INTRODUCTION AND MOTIVATION 

Hadrons containing two heavy quarks (hereinafter: doubly heavy) have the potential to provide deep insights into 
the details of confining forces. Furthermore, these particles offer an opportunity to explore the symmetries which 
emerge in hadrons that contain more than one heavy quark. An interesting possibility is that there is a symmetry 
that relates hadrons with two heavy quarks to hadrons containing a single heavy quark. Savage and Wise were 
the first to propose this, and they used it to relate the hyperfine splittings of baryons containing two heavy quarks 
(doubly heavy baryons, or DHBs) to those of heavy mesons containing a single heavy quark [JJ. We refer to this as 
the superflavor symmetry. Subsequently, this symmetry has been used to relate the hyperfine splittings of tetraquarks 
containing two heavy quarks to those of baryons containing a single heavy quark [2]. 

The fundamental notion of the superflavor symmetry is that, in the infinitely massive heavy quark limit, the 
two heavy quarks in a doubly heavy hadron will bind tightly into a heavy diquark in a color antitriplet, and the 
light degrees of freedom will interact with this stationary, pointlike source of anticolor. As a result the details of 
nonperturbative QCD interactions involving the light degrees of freedom in a QQq system will be identical to the 
corresponding interactions in a Qq system. For the ground state, two heavy quarks of the same flavor will exist in a 
spin triplet due to Fermi statistics. Savage and Wise went on to argue that the interactions between the two heavy 
quarks are dominated by a color Coulomb potential due to the small average separation between them. Under this 
assumption, the energy of the heavy diquark scales as a 2 s (MQ)MQ, in analogy to positronium. This implies that, for 
sufficiently massive quarks, the energy of these excitations will be much larger than the energy required to excite the 
light degrees of freedom. However, several authors have obtained results for doubly heavy baryons which disagree 
with this estimate [3HTTj. Indeed, the Roberts- Pervin (RP) model predicts that, as the heavy quarks become more 
massive, the energy required to excite the heavy diquark decreases. If this is the case, then one cannot ignore these 
excitations, and it is desirable to extend the superflavor symmetry so that it incorporates this possibility. This paper 
seeks to test the validity of such an extension in the context of a constituent quark model We also demonstrate 
the necessity of this extension of the superflavor symmetry by estimating the energy required to excite heavy diquarks 
using experimental data for doubly heavy mesons. 

There are many articles that calculate DHB spectra. Among these treatments arc nonrelativistic constituent 
quark models (NRCQMs) which deal with the three body problem directly [3HS IT2TU9| . There are also several 
NRCQMs which solve the three-body problem through the Hyperspherical formalism [10, 2 23] and some which 
solve the problem through a quark-diquark approximation (we discuss this method in section |II A[ ) [BHEI [23] . There 
are relativistic quark models which treat DHBs, by a quark-diquark approximation [H1[2S][2S] and otherwise [2"Tll30| . 
In addition, DHBs have been treated with bag models (4] [31], the Bethe-Salpeter equation [11] [32] [33] , QCD sum 
rules [T3 1 I34H39 ] . effective field theories [T]l40ti43]. Lattice QCD [44H46]. and by solving Non-Relativistic QCD on the 
lattice [37] [35] . Predictions for the ground states of DHBs have also been made by exploiting regularities in hadron 
spectra and combining those regularities with inequalities derived through the Hcllmann-Feynman Theorem |49l 150) . 

We discuss six aspects of the spectra of doubly heavy hadrons. First, we develop an energy scaling rule for systems 
consisting of two heavy quarks, such as doubly heavy mesons and heavy diquarks, which shows how the excitation 
energy scales with the reduced mass of the system. We then explore a possible relationship between the excitation 
energies of heavy diquarks and those of doubly heavy mesons. The third aspect is the expected existence of degenerate 
multiplets in the spectrum of DHBs, similar to those of Heavy Quark Symmetry. The fourth aspect we discuss is 
the expectation that the energy required to excite the light degrees of freedom of a DHB will be independent of the 
specific excited state of the heavy diquark. The fifth aspect is a heavy diquark flavor symmetry which relates the 



2 



properties of 2 CC , Hh c , and S(,(, baryons. The sixth aspect is a superfiavor symmetry which relates DHBs with a heavy 
diquark in an arbitrary excited state to singly heavy mesons. Finally, we examine how these six aspects emerge (or 
not) in baryons containing strange quarks, with the strange quark treated as a heavy quark. We examine these ideas 
using results from a non-relativistic quark model [5] and, when possible, experimental data. Factorized (see section 



II A), quark-diquark approaches to DHBs should naturally contain these DHB symmetries. However, the model that 



we employ was not constructed with these symmetries and does not assume a quark-diquark structure. Therefore, 
one of the goals of this manuscript is to investigate if these symmetries emerge in such a model. 



II. HEAVY-HEAVY SYSTEMS 



A. Factorization 



We begin our discussion by exploring the physics of doubly heavy mesons and heavy diquarks. In order to accomplish 
this, we work in a limit where an approximation known as factorization [7J 1 1 1 1 15 lj is valid for DHBs. If two quarks 
are sufficiently heavy the potential between them is expected to be essentially Coulombic, and the size of the QQ 
subsystem of a DHB will decrease as Mq increases, with tqq ~ 1/Mq -C 1/A.qcd- Factorization posits that, 
for a pointlikc heavy diquark, one may first solve for the spectrum of the (heavy-heavy) diquark system, and then 
subsequently, solve for the energy of the heavy diquark-light quark (heavy-light) system. Thus, a three body problem 
(in a constituent quark model picture) factorizes into two, independent two-body problems. Unfortunately, the c 
and b quarks may not be massive enough for a diquark composed of these quarks to be considered a pointlike object 
bound by a color Coulomb potential 2 a . However, a pointlike diquark is not a necessary condition for factorization. 
A simple example of this is a system of three particles confined by a harmonic oscillator potential. The Hamiltonian 
factorizes explicitly after a suitable choice of coordinates, provided that the masses of the heavy quarks are equal |52) . 
This approximation complements our discussion about the symmetries of the heavy-light system in that factorization 



can easily incorporate the spin and (super)flavor symmetries that we will present in section III 



Factorization may be viewed as a first approximation to a more careful treatment of the interactions between the 
light degrees of freedom and the heavy diquark. In this approximation, corrections due to the finite size of the heavy 
diquark may be handled using form factors [S] . The validity of this approach also depends on whether confining forces 
are pairwise or three-body; factorization may be expected to be more accurate for three-body rather than pairwise 
confining forces [5J [5T]. For the discussion of factorization, we work in the limit where the ground state and the 
excited states of the heavy diquark arc sufficiently pointlikc so that we may neglect corrections due to its finite size. 

If factorization holds reasonably well for baryons containing b and c quarks, it may be possible to extract the spectra 
of heavy diquarks from experimental data. This may make it possible to investigate the symmetries of a bound state 
of two heavy quarks in an antitriplet color configuration, with the understanding that this color charge is screened 
by the light degrees of freedom. Howbeit, these screening effects do not significantly alter the physics of the heavy 
diquark system if the diquark is sufficiently pointlike. Thus, the interactions between the quarks in the QQ system 
of a DHB should be similar to those between the quark and antiquark of a doubly heavy meson 7\. The primary 
difference arises from the color factors. Naively, these factors are the same for every type of interaction involved, 
implying that the strength of interactions in the doubly heavy meson system are twice that of interactions in the 
heavy diquark system. This is true at tree level in perturbative QCD, and there are indications from quenched latice 
QCD that this is also true of the long-distance part of the confining potential [53]. If this holds rigorously, then it 
amounts to a symmetry between color antitriplet and color singlet bound states containing two heavy quarks. 



B. Interaction Strength and Mass Scaling 



The appropriate effective theory for exploring heavy- heavy systems is nonrelativistic QCD (NRQCD) 40J, but here 
we choose a simpler and more familiar approach. We work in the limit where the two heavy quarks are massive 
enough that their motion can be treated as essentially non-relativistic and QCD interactions can be well-described by 
an adiabatic potential. In other words, we are in the limit where the Schrodinger equation is valid. For simplicity, we 
also ignore all spin-dependent interactions and assume a flavor independent potential. The most general Hamiltonian 
for the excitation energy of two heavy quarks that can be written for this situation is 

H = ~ + bV(r), (1) 

where r is the relative coordinate between the two heavy quarks, ijl is the reduced mass of the system, and P is 
the momentum conjugate to r. We have also included a strength parameter, b, as an explicit factor multiplying the 
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potential, V(r). This Hamiltonian does not mix states with different values of orbital angular momentum because 
there are no scalar spin-independent interactions which can do so. Our goal is not to solve the eigcnenergy problem 
explicitly, but to discern how the eigenenergies scale as one changes the reduced mass and the strength parameter. The 
appropriate tool for accomplishing this task is the Hellmann-Feynman Theorem, which has been applied previously 
to hadron systems (TSJ HH1 EH 58 , but our approach is somewhat novel. 

The Hellmann-Feynman theorem states that, if a Hamiltonian, H, its nth eigenenergy, E n , and the nth eigenstate, 
\^ n ), are differentiable functions of a parameter, A, then the energy satisfies, 

For simplicity, we drop all of the labels referring to the eigenstate and remove all explicit references to the eigenfunc- 
tions. An application of this theorem to the Hamiltonian of eqn. (JlJ leads to 

Si Ul 1 

(V(r)). (4) 



BE 
~db 

On the other hand, the eigenenergy may be written as, 



After a simple substitution, we obtain 



(P 2 ) 

E = ^L + b(V(r)). (5) 



dE BE 

E =-^ + h m- (6) 



Provided that E ^ 0, this partial differential equation is separable. In this case, the solution to the PDE is 

E n i = D nl a d - X b d , (7) 

where D n i is independent of u and 6, n labels the radial excitation, I labels the orbital excitation and d is a constant 
index to be determined, which depends on the specific potential that is used, but is independent of the particular 
eigenstate, labeled by n and I. We show that this is true for a simple case below. This equation was derived previously 
by Quigg and Rosner |54j for the special case of a power law potential. For a potential which is positive definite, 
the solution is exact. However, if the potential, V{ro) — for some Tq € (0, oo), then there may be some value of 
the reduced mass for which the eigenenergy crosses zero. In this respect the solution above can be regarded as an 
approximation where one chooses the value of d so that the equation reproduces the spectrum for a limited range of 
reduced masses. It may be possible to find a more general solution to the PDE above, but this form of the solution, 
eqn. Q, is quite convenient for our purposes. 

C. Scaling for Special Cases 

To illustrate the use of the energy scaling rule (eqn. ny, consider an attractive power law potential, 

V(r) = sgnOK. (8) 
An application of the Quantum Virial Theorem yields, 

(P 2 ) =b»u(V(r)). (9) 
This result can be combined with the results from the Hellmann-Feynman Theorem to obtain, 



The solution to this equation scales as b 2 +» , and so 



d=7^~ (11) 
2 + v K ' 
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for a power law potential. For the special case of a Coulomb potential, the energy is proportional to the reduced mass, 
while the energy for a linear potential scales as and the energy for a harmonic oscillator scales as This 

simple example yields valuable intuition about the meaning of the scaling index, d; if d « | for a range of reduced 
masses then the potential is dominated by a linear term, but if d ~ 2 then the potential is effectively Coulombic. This 
special case for the energy scaling equation of a power law potential agrees with the result which Quigg and Rosner 
obtained using another method [54 . 



D. Doubly Heavy Meson Ratios 



The energy scaling rule, eqn. ([7]), provides a simple way to analyze how orbital and radial excitations of a heavy- 
heavy system scale with the mass of the heavy quarks. Because D n i is independent of the reduced mass, /i, the scaling 
equation predicts that 



M n>l — - M nl - 
QiQ 2 Q1Q2 

M n '}'—, - M n \-, 
Q\Q i 2 Q1Q2 




(12) 



where M nl — is the mass of a doubly heavy meson with radial and orbital quantum numbers n and I, respectively, 
Q1Q2 

and quark flavors Qi an d Q2- MQ1Q2 ^ s the reduced mass of the Q\Q 2 system. This expression is only reliable in 

the nonrelativistic, adiabatic limit, and it contains an inherent ambiguity, in that the reduced mass of the heavy 

quark system cannot be directly measured by experiment. The values of the reduced masses will be model dependent, 

and models which allow for running masses should evaluate each reduced mass at the appropriate renormalization 

scale. Corrections due to the running of the coupling constant may also be warranted. These ratios also ignore 

spin-dependent interactions and, as such, they should be approximately valid for spin-averaged masses. The striking 

prediction that this relation makes is that a host of ratios involving different excited states of doubly heavy mesons 

will be identical: no matter what choices are made for n, n', I, and V the ratio is the same. If radial and orbital 

splittings for doubly heavy mesons satisfy this prediction, then it will support the reasonableness of the energy scaling 

rule, eqn. 



III. HEAVY-LIGHT SYSTEMS 



A. Total Angular Momentum Decoupling 



In this discussion we use the successful Heavy Quark Symmetry (HQS) |59j as a guide for exploring the symmetries 
which might emerge for doubly heavy hadrons. For hadrons containing a single, infinitely massive heavy quark, the 
spin of the heavy quark decouples from the spin of the light degrees of freedom. Thus, singly heavy hadrons exist in 
degenerate doublets with total angular momenta J = J/ ± sq, where sq = \ is the spin of the heavy quark and J/ is 
the total angular momentum of the light degrees of freedom. For doubly heavy hadrons, it is a natural extension of 
this concept to expect that the total angular momentum of the heavy diquark should decouple from the total angular 
momentum of the light degrees of freedom leading to degenerate multiplets with total angular momenta J satisfying 

\Jd- Ji\ < J< Jd + Ju (13) 

where Jd is the total angular momentum of the heavy diquark. If this relation is only applied to situations where 
J^ d — 1 + or + , then one recovers the symmetry proposed by Savage and Wise pQ. However, their formalism may 
easily be extended to diquarks with higher spin |60j . 



B. Flavor and Superflavor Symmetries 



For singly heavy hadrons the excitation spectrum of the light degrees of freedom becomes independent of the heavy 
quark mass, in the infinitely massive heavy quark limit. Thus, the excitation spectrum of a hadron containing a c 
quark is very similar to the excitation spectrum of a hadron containing a b quark. For doubly heavy hadrons, however, 
this symmetry is broken because a consistent effective field theory requires the inclusion of kinetic energy terms for 
the b and c quarks which depend explicitly on their masses [61] . The symmetry that remains for doubly heavy systems 
is the superflavor symmetry l a . If the heavy diquark is indeed a pointlike, stationary source of anticolor, one may 



5 



construct a new effective theory where the relevant degrees of freedom are the heavy diquark and the light degrees 
of freedom. In this theory one may neglect the kinetic energy of the heavy diquark. More specifically interactions 
which change the velocity of the heavy diquark are suppressed by its mass, M d ^> Aq C d- 

In analogy with the heavy quark effective theory (HQET), it is expected that the effective Lagrangian for the 
quark-diquark system will be independent of the diquark's mass. Thus, if a diquark of mass Md and total angular 

momentum and parity J£ d is replaced with a diquark of mass M' d and total angular momentum and parity j7 d , the 
physics of the light degrees of freedom will be unchanged. Additionally, the fact that the effective Lagrangian does not 
depend on the heavy diquark's spin (also in analogy with HQET) indicates that one may also replace a heavy diquark 
(of any angular momentum and parity) with a heavy antiquark, and the physics of the heavy-light system will be 
unaffected. In practice this is a superflavor symmetry which relates the spectra and properties of singly heavy mesons 
to those of DHBs, independent of the excitation state of the diquark. This is the generalization of the superflavor 
symmetry to which we alluded earlier. 

Given the arguments above, a new perspective on HQS appears: Heavy Quark Symmetry may be a special case 
of a larger symmetry between hadrons containing Nq heavy quarks and/or antiquarks in a color configuration C, 
interacting with Ni light quarks in a color configuration C, and hadrons containing Nq heavy quarks and/or antiquarks 

in a color configuration C, interacting with Ni light quarks in a color configuration C. In short, we speculate that 
the decoupling of the spin and flavor/mass of the heavy degrees of freedom from those of the light degrees of freedom 
may be a general feature of QCD in the infinitely massive heavy quark limit. The validity of this generalization of 
HQS can be tested with heavy exotic hadrons, if they exist [2]. 



IV. SYMMETRY CONSTRAINTS 



A. Meson-Baryon Ratios 



If there is a symmetry between color singlet and color antitriplet bound states, then the spectra of DHBs and the 
spectra of doubly heavy mesons can be related in a simple way. Factorization of DHBs implies that one may write the 
mass as the sum of the energy contributions from the heavy diquark and the light degrees of freedom, M = Ed + Ei, 
and after combining this relation with the energy scaling rule eqn. Q, we find that 



— M nl — 
Q1Q2 Q1Q2 

71 jn'l'vX n/fnlvX 

m QiQ2i m QiQ 2 q 



(14) 



where M ™ — is the mass of a doubly heavy meson with radial and orbital quantum numbers n and I, respectively, 

Q1Q2 

and Mq"q is the mass of a DHB which contains a heavy diquark in a state labeled by n and I, and light degrees of 
freedom in a state labeled by v and A. The eigenstate of the light degrees of freedom for the "primed" DHB must be 
the same as that of the "un-primed" DHB in order to cancel its contribution to their masses in the denominator of eqn. 
( |14| . This result assumes that the strength of the interaction between the heavy quark and heavy antiquark in the 
doubly heavy meson is exactly twice the strength of the interaction between the two heavy quarks in the DHB (Sec. 



II A). We have neglected the running of the coupling constant as well, but interactions in the heavy diquark system 
and interactions in the doubly heavy meson system should occur at a similar renormalization scale. The expression 
also neglects spin interactions. We also remark that, inasmuch as factorization is valid, this relation may be applied 
to Qqqi as well as Sqq/ baryons. Because this constraint allows one to compare properties of a meson to a baryon 
through their color structures, the combination of factorization with the symmetry between singlet and antitriplet 
color configurations might be thought of as a "supercolor" symmetry. 



B. DHB Spectrum 



The extension of the superflavor symmetry from section MB stipulates that the physics of the light degrees of 
freedom for a DHB is independent of the flavor and excited state of the heavy diquark. This means that splittings 
between the ground state and an excited state of a singly heavy meson should be the same as the analogous splitting 
in the light degrees of freedom of a DHB. The DHB spectrum should consist of a pattern of states that repeats for 
each heavy diquark with total angular momentum and parity J^ d . Furthermore, this repeating pattern will be nearly 
identical to the singly heavy meson spectrum, after corrections due to spin interactions and the finite heavy quark 
mass have been made. The decoupling of the total angular momentum of the heavy diquark will also lead to a DHB 
spectrum consisting of degenerate multiplets with total angular momentum, | J<j — J{\, + J/. 



G 



FIG. 1: DHB Symmetries. A DHB with two heavy quarks, Q and a light quark q can be treated as two independent systems 
due to factorization and the decoupling of the total angular momentum of the heavy diquark d. The spectrum of the heavy-light 
d — q system can be related to the spectrum of heavy mesons through the superflavor symmetry. Additionally, the d — q system 
can be related to a d! — q (where d! is a heavy diquark of a different flavor) system through the diquark flavor symmetry. The 
spectrum of the heavy-heavy Q — Q system can be related to the spectrum of doubly heavy mesons through the supercolor 
symmetry. 



\ 

V 



q) 

y 



Factorization 



Diquark Total Angular 
Momentum 
Decoupling 



< 



d 




q 

d' 



fn \ Supercolor Symmetry \ 

\ oh q) 

v y v y 



V. QUARK MODEL FOR DHBS 



A. Hamiltonian 



In this section we present a non-relativistic quark model which Roberts and Pervin [3] used to calculate masses 
of heavy baryons. The advantage of this model is that it is among only a handful of articles [H [3 [3 HH] m which 
orbital excitations both in the heavy diquark and in the light degrees of freedom are allowed. Among these models, 
it is the only one which solves the three-body problem without using a quark-diquark approximation. The model is 
also convenient for this manuscript because we have access to the wave functions and can comment on the multiplet 
structure of the states obtained in the spectrum. For the purposes of this note, the spectrum of [3] has been recalculated 
with spin-dependent interactions removed. The Hamiltonian is 

H = K + V con f + C qqq . (15) 

C qqq is a constant term which is the same for all baryons in the model. The kinetic energy term is 

2 



i=l,3 v l/ 

and it may be factorized easily through the introduction of Jacobi coordinates, 

1 



and 



^(ri-r 2 ) (17) 



A= /2/ m iri+ m 2 r 2 _ r (lg) 



3 V mi + ni2 
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With this choice, p is proportional to the relative coordinate between the two heavy quarks, and A is proportional to 
the relative coordinate between the center of mass of the heavy diquark and the the light quark. In the rest frame of 
the baryon, the kinetic energy operator becomes 



K=y i m i + ^ + ^, (19) 



Pi , Px 



=1.3 



with /i p = ^"^ 2 2 and p\ — f (|™ 1 H ^^)1^ 3 ■ This particular model contains pairwise interactions for the confining 
potential, 

V coni = ± (*f* - ^ , (20) 



»<?'=! 
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with Tij — \ri — Yj\. The ri2 terms combined with the p term in the kinetic energy give a Hamiltonian for the 
heavy- heavy subsystem which is in the form of eqn. 0, the form we used earlier to derive the energy scaling rule, 
eqn. (f7|. Unfortunately, the terms in ri 3 and r 2 3 depend on both p and A. The model therefore does not factorize 
explicit y, nor is there any obvious limit in which factorization emerges. 

Obtaining a good fit to the baryon spectrum for this model required setting the strength of the Coulomb parameter, 

acoulj to be quite small (~ 1CT 4 ). As such, we expect the eigenenergy contribution from the "p" part of the 

_ i 

Hamiltonian to scale roughly as p p 3 , which explains why the model predicts that the energy required to excite the 
heavy diquark decreases as the heavy quark mass increases. The small value of acoui also has implications for the size 
of heavy diquarks in this model. For a Coulomb-plus-linear potential one expects that the size of the heavy diquark 
approaches tqq ~ 1/Mq as the quarks become very heavy since the Coulomb term dominates. Similarly, the linear 
term dominates as the two quarks that comprise the diquark become very light, and the diquark size approaches 
tqq ~ l/(AfgaLin) 3 ■ Because the potential in this model is essentially linear, the RP model might be considered 
a worst-case scenario for factorization because it contains pairwise forces and because the heavy diquarks are less 
pointlike than those of other Coulomb-plus-linear NRCQMs. We have listed the parameters of the model in table [T] 



TABLE I: Parameters of the RP model j3j . m a is the mass of the u and d quark. 



m a m 3 m c m b a^in acoul C qqq 
(GeV) (GeV) (GeV) (GeV) (GeV 2 ) - (GeV) 
0.2848 0.5553 1.8182 5.2019 0.1540 « 0.0 -1.4204 



B. Flavor Multiplets 

Assuming that only the b and c quarks can be treated as heavy, states containing two heavy quarks can be placed 
in SU(2) multiplets. These multiplets are a triplet 

E cc = ccq, E bb = bbq, E bc = (cb + be) q, (21) 



and a singlet 



E' bc =^=(cb-bc)q, (22) 



with the light quark q being one of u, d or s. Of course, this is only meant to be a classification scheme, but it has 
implications for the angular momentum structure of the states, since the spin-space-flavor wave function must be 
symmetric in 'identical' quarks. For instance, if there is no excitation of the heavy diquark, the diquarks in the triplet 
must have spin 1, while those in the singlet must have spin 0. 



C. Angular Momentum Supermultiplets 

In the quark model presented above, states arc defined in what we may call an L-S basis. Symbolically, the angular 
momenta in the states are coupled in the scheme {[l p l\]L[sdSi}s}j, where l p is the relative orbital angular momentum 
between the two heavy quarks, l\ is the relative orbital angular momentum between the light quark and the heavy 
diquark, Sd is the spin of the heavy diquark, and s; = ^ is the spin of the light quark. For the symmetry scheme 
discussed in Section 



III 



states are identified by the spin and parity of the diquark, J d d , the spin and parity of the 
light degrees of freedom, J,* 7 , and the total angular momentum and parity of the baryon, J p . In order to compare 
these states with quark model states it is necessary to transform the states from the L-S basis to the Jd~Ji basis, 
in which the angular momenta are coupled symbolically as {[^ p Sd] j d [l\si] j, }j. These two bases are related by a 9-j 
symbol, 



lp 


h 


L 




1 

2 


S 


Jd 


Jl 


J 



V(2£ + 1)(25 + l)(2J d + 1)(2J, + 1) ^ s d i S }. (23) 



Using this, we can assign states obtained in the quark model to angular momentum supermultiplets by examining 
their wave function components. 



D. Pair- wise Mixing 

The extension of the superflavor symmetry that we have proposed is broken in this model even when spin-dependent 
interactions are absent. This arises because of pair- wise forces and the composite nature of the heavy diquark, and it 
is a feature of any model that includes such forces. To illustrate, let us consider a scalar, pair-wise interaction /(|rij|) 
between quarks i and j. Such an interaction may be written quite generally as, 

oo 

f(\nj\) = X>(p, A)(y,(ty,) • r,(fijO). (24) 

1=0 

This expression comes from writing ffa and r*23 as linear combinations of p and A and expanding the function / as 
a sum over spherical harmonics multiplied by a scalar function gi. For a Coulombic interaction, this equation yields 
the familiar multipole expansion. Matrix elements of the interaction / are proportional to a 6 — j symbol, 



(I'l^Lim-rMphlL)*! »?. (25) 




This operator, eqn. (24), is a scalar in the total orbital angular momentum of the system, L. It therefore mixes 
states having different values of l p and l\, but the same total L. For instance, states having l p = 0, l\ = 2 mix 
with states having l p — 2, l\ = 0. However, even if we ignore mixing between states with different values of l p and 
l\, the interaction is diagonal in the L-S basis, implying that it mixes states in the Jd~Ji basis. In the absence of 
other interactions, the mixing coefficients are independent of the mass of the heavy quarks and the strength of the 



interaction, and are given by eqn. (23). In general, this mixing also leads to a mass splitting between states that 
occupy the same symmetry multiplet, and thereby lifts the degeneracy of states within the multiplct. However, if 
l p = I' = or l\ = l' x = 0, the matrix elements of pairwise operators simplify because the 6-J symbol of eqn. (25) 



dictates that I — 0. The matrix element of the interaction is then diagonal in the Jd~Ji basis, as well as the L-S 
basis. The upshot is that these pairwise interactions only remove the degeneracies within multiplets which contain 
both an orbitally excited heavy diquark and orbitally excited light degrees of freedom. 



VI. RESULTS 



It must be emphasized that the majority of treatments of DHBs assume factorization, so that most of the features 
that we have discussed arise naturally. For instance, the fact that the excitation energies of the light quark are inde- 
pendent of which state the heavy diquark is in, is a natural consequence of factorization. So, too, is the independence 
of these excitation energies of the flavor of the diquark. The angular momentum 'supermultiplets' should also arise 
naturally in factorized models. Our goal is to elucidate whether or not models that do not assume factorization, or, 
to be more precise, how a particular model that does not assume factorization, manifests these predictions. 
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A. Meson-Meson Ratios 



The meson-meson ratios of eqn. (12 1 provide an opportunity to determine the energy scaling index d using experi- 
mental results. There are almost enough experimental data to compute the spin averaged masses of the ground state, 
the first orbital excitation, and the first radial excitation for charmonium and bottomonium |62j . For the hb{lP), we 
use the recent result from the Belle Collaboration [53]: for the ground state bottomonium, we use the measurement 
from the BaBar Collaboration [64]; and for the r)j,(2S), we use the result from the model of Ebert et al. [65] • With 
these data there are two ratios that can be constructed: the ratio of the orbital splitting of charmonium to that of 
bottomonium, and a corresponding ratio for radial splittings. The results are shown in tabic [TTJ The energy scaling 
rule, eqn. ([7]), predicts that these two ratios should be identical, and they are up to 6%. This suggests that the 
assumptions behind the derivation of this scaling relation (non-rclativistic dynamics, constant scaling index for a 
finite range of reduced masses), are reasonable. The experimental ratios also encapsulate a non-trivial fact about 
the spectra of doubly heavy mesons - the cost of an orbital or a radial excitation is nearly independent of the mass 
of the two heavy quarks. This fact forces non-relativistic models to choose a potential such that the energy scaling 
parameter, d 1 for the reduced masses of charmonium and bottomonium. This is quite different from the energy 
scaling for a Coulomb potential, and contradicts some formulations of NRQCD which use a power counting scheme 
that assumes Coulombic excitation energies [35] . Quigg and Rosner have shown previously that the energy splittings 
of doubly heavy mesons can be reproduced using a suitably tuned Coulomb-plus-linear potential or a logarithmic 
potential [54] . 



TABLE II: Spin-averaged mass splittings for bottomonium and charmonium from experiment [62] • The energy scaling rule, 
eqn. (7 1, predicts that the ratios of these splittings should be the same value, ^^"^ • This allows us to extract the energy 
scaling index, d, from these ratios using the heavy quark masses from the RP model [3]. 



Splitting Am fi (MeV) Am s5 (MeV) Amg/Am^ d 
M(1P) - M(1S) 457.6 ± 0.31 457.0 ± 1.5 1.001 ± 0.003 0.999 ± 0.007 
M(2S) - M (IS) 606.1 ± 1.0 573.3 ± 1.6 1.057 ± 0.003 0.947 ± 0.007 



B. Meson-Baryon Ratios 



The meson-baryon ratios from eqn. ( 14 1 provide a means of comparing DHBs and doubly heavy mesons, and these 
ratios are sensitive to the form of the confining potential. We use this ratio to compare two different "models" . We 
construct ratios involving DHB masses from the RP model [3] and the results of a simple model for doubly heavy 
mesons using the Hamiltonian, eqn. |l]), with its parameters and potential set to match the RP model (see table 



Illl. The RP model violates the meson-baryon ratios, eqn. ( [14] ), for a linear potential by 23% or less. However, their 

potential is not purely linear in the heavy diquark's relative coordinate, and the values for the ratios in the table 

above are all near 1.26. This suggests that effects due to the three-body nature of the DHB system have modified the 

effective interaction strength of the QQ subsystem by a factor of 1.41. This may be due to screening effects from the 

light quark, or it may be a geometrical factor arising because of the pairwise confining potential used in the model. 

The meson-baryon ratios also provide a convenient way to estimate the orbital and radial splittings of heavy 

diquarks. Using d = 1 for the scaling index, we find that the radial splitting for heavy diquarks should be roughly 

300 MeV and the orbital splitting should be about 230 MeV. Because the energy scaling index d « 1 these estimates 

do not depend on quark model masses, and in this sense they are model independent. These estimates assume that 

screening effects from the light degrees of freedom are negligible. For the RP model the radial splittings of the S cc , 

Etc, and Ebb, diquarks are 332, 294, and 223 MeV, respectively; likewise, the orbital splittings are 188, 167, and 129 

_ i_ 

MeV. This disagreement with the doubly heavy meson spectra is due to the \i d 3 scaling of excitation energies in the 
RP model. 
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TABLE III: Meson-Baryon ratios (eqn. 



14 1 from the spin independent RP model and a simple quark model for doubly heavy 

The subscript m refers to a meson and the 
23 = 1.59. 



mesons. For all of these ratios the light degrees of freedom are in the ground state 
subscript, b, refers to a baryon. The ratios of these splittings should be the same value 



Ratio (ridld) cc be bb 



M m {lP) - M m (lS) 
M 6 (1P) - M b (lS) 

M m (2S) - M m (lS) 
M b (2S) - M b (lS) 

M m (lD) - M m (lS) 
M b (lD) - M b {lS) 



1.28 1.26 1.32 
1.24 1.23 1.30 
1.27 1.25 1.33 



C. Symmetry Multiplets 

The quark model states for the 2 CC , H& c , and the from the spinless RP model are arranged into symmetry 
multiplets in table |IV| and the results from the full model [5] are in table |V| The results for the full model and the 
spinless model for flavor-singlet E' bc states are in table 

interactions lift the degeneracy of the states by less than Yt> MeV, and these splittings decrease as the heavy diquark's 
mass increases. One interesting feature of the model is that there are degeneracies in the spectrum in addition to those 
predicted by the decoupling of the heavy diquark's spin. These additional degeneracies arise due to the fact that the 
RP model does not include a fully microscopic spin-orbit interaction. For instance, in table|v]the n 2 , Sd+1 (l p )j d = 1 3 -Di, 
1 3 I?2, and 1 3 Z?3 multiplets are nearly degenerate. An appropriate spin-orbit interaction would split these states, but 
this would be suppressed by the masses of the heavy quarks. There are also additional degeneracies among the 
n 2 d Sd+1 (l p ) Jd , n x {l\) Jl = l 3 Si, lPi and l 3 ^, IPs multiplets as well as the 1 1 P 1 , lPi and l l P u IPs multiplets. 
These degeneracies would be lifted if the model included a spin-orbit interaction which couples the orbital angular 
momentum of the heavy-light system to the spin of the light quark, such as an operator proportional to l\ • s/. Such 
an interaction is not expected to be suppressed in the heavy quark limit. 



VI 



Tables [V] and | VI| show that the inclusion of spin dependent 



TABLE IV: Symmetry multiplets for the RP model for flavor-triplet states with spin interactions removed. The decoupling of 
the diquark total angular momentum predicts that the states within a multiplet will be degenerate. 



Multiplet 


J 


P 






Mass (GeV) 




n d Sd+1 ( l p)jd> n ^i l ^Ji 


(\Jd-Ji\,- 


■ ■ ^ 


Jd + Jl) 




Hbc 




l 3 Si, lSi 

2 




3 
2 


f ) 


(3.724, 3.724) 


(7.045, 7.045) 


(10.341, 10.341) 


I 1 Pi, ISi 

2 




3 
2 


") 


(3.912, 3.912) 


(7.212, 7.212) 


(10.470, 10.470) 


2 3 S ly ISi 

2 




3 
2 


f ) 


(4.056, 4.056) 


(7.339, 7.339) 


(10.564, 10.564) 


1 3 D U ISi 


( i + 
V2 ' 


3 
2 


f ) 


(4.078, 4.078°) 


(7.358, 7.358°) 


(10.579, 10.579°) 


1 3 D 2 , ISi 


( 3 + 


5 
2 


f ) 


(4.078°, 4.079) 


(7.358°, 7.360) 


(10.579°, 10.581) 


1 3 D S , ISi 

2 


(f + , 


7 
2 


f ) 


(4.079, 4.078) 


(7.360, 7.358) 


(10.582, 10.579) 


l 3 Si, lPi 

2 


(!"- 


3 
2 


~) 


(4.073, 4.073 6 ) 


(7.391, 7.391°) 


(10.691, 10.691°) 


l 3 5i, IPs 


(I" 3 
V2 ' 2 


) 


I") 


(4.073, 4.073\ 4.073) 


(7.391, 7.391°, 7.391) 


(10.691, 10.691°, 10.691) 


r^Pi, ipi 


( i + 
V2 ' 


3 
2 


f ) 


(4.259 c , 4.249) 


(7.557 c , 7.548) 


(10.818 c , 10.813) 


V-Pl, IPs 

2 


(1+ 3 
^2 ' 2 


+- 

) 


f + ) 


(4.235 c , 4.259, 4.247) 


(7.536 c , 7.557, 7.546) 


(10.804 c , 10.818, 10.809) 


l 3 Si, 2Si 

2 


f i + 

V2 ' 


3 
2 


f ) 


(4.346, 4.346) 


(7.662, 7.662) 


(10.964, 10.963) 


l 3 ^!, IDs 

2 


(1+ 3 
V2 ' 2 


4- 

) 


2 > 


(4.392, 4.392, 4.385 d ) 


(7.709, 7.709, 7.700 d ) 


(11.013, 11.013, 11.001 d ) 


l 3 ^!, ID 5 

2 


(3+ 5 
^2 ' 2 


+- 

) 


2 > 


(4.392, 4.385 d , 4.392) 


(7.709, 7.700 d , 7.709) 


(11.013, 11.001 d , 11.013) 



abcd These states mix significantly. 
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The energy splittings due to pairwise mixing (section VD) are small in the RP model. For instance, the mass 
splittings due to the mixing of states with l p = 0, l\ — 2 and l p — 2, l\ — are not noticeable in the spectra. The 
mixing of states due to this is also quite small and vanishes as the heavy quark mass increases (although there is no 
a priori reason to expect it to). States that are orbitally excited in both the heavy diquark and the light degrees of 
freedom are also subject to pairwise mixing in the Jd and J/ quan tum numbers. This is the reason for the broken 
degeneracies of the n d Sd+ (l p ), 
column of table 



VI 



substantial, particu 



n\{l\) = 1 1 P, IP multiplets in table IV and the 1 3 P, IP multiplcts in the "Spinless" 



All of these splittings are small (less than 25 MeV) , but the mixing associated with them is 
arly for the flavor-singlet states. 



TABLE V: Symmetry multiplets for the RP model for flavor-triplet states including spin dependent interactions. The decoupling 
of the diquark total angular momentum predicts that the states within a multiplet will be approximately degenerate. 



Multiplet 

nf d+1 (l P )j d , n x (l x )j l 


J p 

(| Jd — Ji\, Jd - 


- Ji) s cc 


Mass (GeV) 


Hf,6 


l 3 5i, ISi 

2 


\2 ' 2 > 


(3.678, 3.752) 


(7.014, 7.064) 


(10.322, 10.352) 


I 1 Pi, ISi 

2 


,1- 3-N 

y-2 ' 2 > 


(3.911, 3.917) 


(7.212, 7.214) 


(10.470, 10.470) 


2 3 5i, ISi 

2 


n+ 3+s 

V2 ' 2 > 


(4.030, 4.078) 


(7.321, 7.353) 


(10.551, 10.574) 


l 3 Di, ISi 


V2 ' 2 I 


(4.098, 4.045 a ) 


(7.372, 7.337°) 


(10.589, 10.564°) 


1 3 D 2 , ISi 

2 


/3 + 5+\ 
V2 ' 2 > 


(4.094°, 4.092) 


(7.369°, 7.368) 


(10.587°, 10.588) 


i 3 £> 3 , ISi 

2 


V2 ' 2 > 


(4.048, 4.095) 


(7.340, 7.370) 


(10.568, 10.588) 


1 3 S U lPi 


/I- 3-N 
V2 ' 2 I 


(4.081, 4.077) 


(7.397, 7.392) 


(10.694, 10.691) 


l 3 Si, IPs 


(1- 3" S-S 
^2 '2 '2 > 


(4.073, 4.079, 4.089) 


(7.390, 7.394, 7.399) 


(10.691, 10.692, 10.695) 


1'PI, lPi 

2 


fl+ 3-K 
V2 ' 2 1 


(4.257\ 4.253) 


(7.555°, 7.549) 


(10.817\ 10.812) 


PPl, IPs 

2 


(1+ 3+ 5+\ 
V 2 '2 '2 / 


(4.230°, 4.261, 4.259) (7.534 6 , 7.557, 7.549) (10.802 6 , 10.818, 10.811) 


l 3 Si, 2Si 

2 


\2 ' 2 


(4.311, 4.368) 


(7.634, 7.676) 


(10.940, 10.972) 


I 3 Si, 1P>3 


/1+ 3+ 5+\ 
^2 '2 '2 J 


(4.394, 4.394, 4.387) 


(7.709, 7.708, 7.701) 


(11.011, 11.011, 11.002) 


I 3 Si, l-Ds 

2 


/3+ 5+ 7+\ 
^2 '2 '2 ' 


(4.391, 4.388, 4.393) 


(7.706, 7.702, 7.708) 


(11.011, 11.002, 11.011) 



These states mix significantly. 



D. Diquark Excitation Symmetry 



The extended superflavor symmetry (Section IIIB I, dictates that the DHB spectrum will exhibit a repeating pattern 



corresponding to the excitations of the light degrees of freedom; each occurrence of this pattern will be "shifted" by 



the mass of the heavy diquark. Table VII shows the energy required to excite the light degrees of freedom in flavor- 
triplet DHBs in the spinless RP mod el for multiplets with two different diquark excitations: the ground state and the 



2s d + l, 



(l p )j d =l 1 Pi multiplets. Table VIII shows analogous splittings for flavor-singlet states, f or th e ground state and 



the Pj d diquarks. The excitation symmetry requires that each entry in the first row of table VII should be equal to 



the entry beneath it on the second row, and each entry on the third row should be identical to the one beneath it on 



the fourth row. Likewise, the entries in the first four rows of table VIII should be equal, as should the entries in rows 
five to eight. This symmetry should be broken by this particular model because of pair-wise interactions, but all of 
the excitation energies that should be identical agree to within 14 McV for the flavor-triplet states, and 16 MeV for 
the flavor-singlet states. We can therefore conclude that the effects of pair-wise forces on this aspect of the spectrum 
are small. 
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TABLE VI: Symmetry multiplets for the RP model for flavor-singlet states. The decoupling of the diquark total angular 
momentum predicts that the states within a multiplet will be approximately degenerate. 



Multiplet 


J p 




Mass (GeV) 


n 2 d d+1 {lp)j d , n x (l\).j l 


(\Jd — Jl\, Jd + Jl) 


Spinless 


Spin 


l x S , lSi_ 


(| + ) 


(7.053) 


(7.037) 


1 3 P , ISi 

2 


(D 


(7.219°) 


(7.230°) 


l 3 Pi, 15 1 

2 


(1- 3-N 
W ' 2 / 


(7.219 a , 7.219) 


(7.199°, 7.228) 


1 3 P 2 , 15i 


(3- 5-N 
V2 ' 2 I 


(7.219, 7.217) 


(7.201, 7.265) 


2 1 S , ISi 




(7.349) 


(7.333) 


l 1 D 2 , ISi 


(IM + ) 


(7.367, 7.367) 


(7.361, 7.367) 


l x S Q , IPi 


(D 


(7.398) 


(7.388) 


I 1 So, IPs 

2 


(I") 


(7.398) 


(7.390) 


l 3 Po, IPi 

2 




(7.553°) 


(7.548°) 


l 3 Po, IPs 

2 


(l + ) 


(7.561) 


(7.546°) 


l 3 Pi, IPi 

2 


V 2 '2 J 


(7.562\ 7.553 c ) 


(7.552\ 7.556°) 


I 3 Pi, IPs 


(1 + 3+ 5+\ 
\2 '2 '2 J 


(7.562 6 , 7.553 c , 7.553 d 


) (7.555 6 , 7.548°, 7.553 d ) 


1 3 P 2 , IPi 

2 


(•3+ 5+\ 
\2 ' 2 J 


(7.541 c , 7.553 d ) 


(7.563°, 7.548 d ) 


1 3 P 2 , IP3 

2 


(1+ 3+ 5+ 7+\ 
V 2 '2 '2 '2 / 


(7.541 6 , 7.561 c , 7.561 d , 7.553) (7.519°, 7.562 c , 7.559 d , 7.557) 


l^o, 25i 

2 


(| + ) 


(7.661) 


(7.645) 


l^o, IDs 

2 


(| + ) 


(7.707) 


(7.709) 


l^o, ID 5 

2 


(| + ) 


(7.707) 


(7.689) 



abcd These states mix significantly. 



E. Diquark Flavor Symmetry 



The heavy diquark flavor symmetry dictates that the energy required to excite the light degrees of freedom should 
not depend on the mass of the heavy diquark. This can be seen for flavor-triplet states in tabic \VU\ where the entries 
on each row are approximately equal: in the spinless Roberts-Pervin model these relations work out within 7 MeV. 
We wish to re-emphasize that this result has emerged from a model which does not explicitly factorize in the heavy 
quark limit and does not assume anything about the s ize of the heavy diquark. We can also see if the e nerg y splittings 
depend on the statistics of the heavy diquark. Table VIII shows the same energy splittings as table VII 
states. As an example, the energy required to excite the light degrees of freedom from the n\(l\)j l = IS 

(rows 1 and 2) should match the excitation energies in table 
states with the splittings for the S& c states we see that 



for the S' 



he- 



wn 



to the IPi multiplets in table 
4). Comparing the splittings for the rl b 
holds within 21 MeV. 



multiplet 
(rows 1 through 



VIII 

the expected relation 



F. Superflavor Symmetry 

In this section we use the superflavor symmetry to compare the spectra of singly heavy mesons to DHBs pQ . We 
do not extensively compare these two types of hadrons here, but rather, we simply compare the energy required to 
orbitally excite the light degrees of freedom of a DHB to the analogous splitting in single heavy mesons. There are 
not enough experimental data to compare with Jp — h states, but there are experimental data for the Jp = | 
states. The flavor symmetry of HQS predicts that these splittings should be the same for D mesons and B mesons 
(neglecting 1/Mq corrections), and the superflavor symmetry dictates that this excitation energy should be identical 



to the third and fourth rows of splittings for DHBs in table VII Table IX contains these splittings for heavy mesons. 
This treatment seems to show that the model of Ref. [3] underestimates the size of the orbital splitting of the light 
degrees of freedom. We also see an indication that the charm quark's mass may still be too light for the flavor 
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TABLE VII: Energy required to excite the light degrees of freedom of a flavor-triplet DHB from the ground state, nx(l\)j, w i = 
+ , to the excited state listed in the second column, for different states of the heavy diquark, listed in the first column. We 

2 

have averaged over the masses in each multiplet. This excitation energy should be independent of the spin, parity, and flavor 
of the heavy diquark, so that the energies in the first row should be the same as the corresponding ones in the second row, and 
the entries in the third row should equal those in the fourth row. These excitation energies should also be independent of the 
flavor of the quarks that comprise the heavy diquark, so that all entries on a particular row should be equal. This is seen not 
only for the first four rows, but also for the last three rows. 



Diquark state 




Upper Multiplet 


Lower Multiplet 


m* — 


m (MeV) 


n d Sd+1 (l P )j d 


n\(h)j l 


(m* 


) 


(m ) 


mcc 


Sic 






lPi 

2 


(§"> 


5 ~) 
2 I 


fl + 3+-, 
V2 ' 2 > 


349 


346 


350 


l X Pl 


lPi 

2 


+ , 


3+ } 
2 I 


/l- 3-N 
V2 ' 2 / 


342 


341 


346 


1 3 S 1 


IPs 

2 


n- s- 

V2 ' 2 


2 > 


/1+ 3+1 
V2 ' 2 > 


349 


346 


350 


l X Pl 


IPs 

2 


(1+ 3 + 
V2 ' 2 


5+) 
2 I 


V2 ' 2 J 


335 


334 


340 


l 3 £l 


ID 3 

2 


fl + 3 + 
<-2 ' 2 


5+) 
2 > 


(1+ 3+1 
\2 ' 2 ' 


666 


661 


668 




ID 6 

2 


(3+ 5 + 
V2 ' 2 


2 > 


f l+ 3+N 
V2 ' 2 > 


665 


660 


666 


l 3 £l 


25*1 

2 


( 1 + 


2 > 


fl + 3+1 
^2 ' 2 J 


622 


617 


623 



TABLE VIII: Energy required to excite the light degrees of freedom of the flavor singlet DHB from the ground state, n\(l\)j^i = 
1S'i + , to the excited state listed in the second column, for different states of the heavy diquark, listed in the first column. We 

2 

have averaged over the masses in each multiplet. The excitation energies on the first four rows should be equal, as should those 
on rows five to eight. 



Diquark state 
(l P )j d 



2s d + l 



l 3 Po 

l 3 Pl 
1 3 P 2 

l'So 
l 3 Po 

1 3 P1 
1 3 P 2 



lPi 

2 

1-Pi 

2 

1-Pi 

2 

lPi 

2 

IPs 

2 

IPs 

2 

IPs 

2 

1P 3 



Upper Multiplet 

(D 



V 2 
( 



3+1 
' 2 ^ 
5+x 



3 

2 ' 2 

(| + ) 
H 3+ 5+\ 
k2 '2 '2 / 
cl+ 3+ 5+ 7+\ 
v2 '2 '2 '2 - 1 



fl- 
^ 2 



Lower Multiplet 

(m ) 

(| + ) 

(D 

fl- 3-1 
V2 ' 2 / 
f3" 5-1 

V 2 '2 / 

(| + ) 

(D 

/l- 3-1 

V 2 '2 



f2" 
1 ■> 



5" 
' 2 



m (MeV) 

5j, c 

345 
334 
339 
329 
345 
342 
337 
336 



l'So 
l 1 ^ 
l'So 



1P>3 

2 

1P>5 

2 

2Si 



(| + ) 
(I + ) 
(| + ) 



(| + ) 
(| + ) 

(D 



654 
654 
608 



symmetry to emerge to a high degree of accuracy. Nevertheless, we can use these results to estimate this splitting to 
be on the order of 400 MeV for DHBs. 

In the RP model there is a separation in mass between states that involve excitations of the heavy diquark and 
states that involve excitations of the light degrees of freedom. For instance, the difference in energy between the 
ground state and the first orbital excitation of the heavy diquark is about 130 MeV for the "Ebb- In contrast, the 
difference in energy between the ground state and the first orbital excitation of the light quark is about 350 MeV. 
Furthermore, the mixing between states with different orbital excitations of the diquark is small. We can determine 



if this separation holds using the estimate of 230 MeV for the orbital splitting of heavy diquarks from section VI B 



and the estimate of 400 MeV for the energy required to excite the light degrees of freedom that we presented above. 
This establishes a hierarchy of energy scales: A d <§; Hq 1 q 2 , A; <C M d ~ Mq 1 + Mq 2 , and A d < A ; ; A d is the energy 
required to excite the heavy diquark, A; ~ Aqcd is the energy required to excite the light degrees of freedom, M d is 
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the mass of the heavy diquark, and ^q 1 q 2 is the reduced mass of the heavy diquark subsystem. 



TABLE IX: Spin averaged masses and the energy (MeV) required to excite the light degrees of freedom of singly heavy mesons 
from the ground state Jj l — | + state to the Jj 1 = | state [62]. We have also included results for Kaons. 



Flavor 


"1(0", I") 


m(i+,2+) 


Am 


K 


793.15 ± 0.20 


1368.0 ± 2.5 


574.9 ± 2.5 


D 


1971.43 ± 0.13 


2446.46 ± 0.52 


475.04 ± 0.53 


B 


5313.70 ± 0.40 


5735.7 ± 2.9 


422.0 ± 2.9 



G. Strange Diquarks 

In many quark models, the strange quark's mass is similar in magnitude to Aqcd- This should be much too light 
for this quark to be considered heavy, but features in the calculated spectrum of ref. [3] emerge that are very easily 
understood if this odd approximation is made. If the strange quark is treated as heavy, the flavor multiplets of section 
VB need to be modified. Instead of the broken SU(2) of that section, we now consider a (badly) broken SU(3) heavy 
Savor symmetry. There will now be a sextet of symmetric states 

3 = ssq, E cc = ccq, E bb = bbq, E c = (sc + cs) q, E b = i (sb + bs) q, E bc = (cb + be) q, (26) 

and an antitriplet of antisymmetric states 

s c = (sc - cs) q, E' b = i (sb - bs) q, E' bc = (cb - be) q. (27) 

For these states, there are only two choices for the light quark, q, namely u and d. These multiplets should not be 
confused with the multiplets in which singly-heavy baryons are placed, with the strange quark taken as one of the 
light triplet of quarks. The SU(3) light sextet then contains the states 

Sq = uuQ, ddQ, —= (ud + du) Q; E'q — —= (su + us) Q, —= (sd + ds) Q; D,q — ssQ, (28) 
v2 V2 v2 

and the corresponding antitriplet consists of 

Aq = (ud - du) Q, E Q = -^=(us- su)Q, -^=(ds - sd))Q. (29) 

In these multiplets, the heavy quark Q may be either b or c. 

In the following subsections we examine the consequences of treating the strange quark as a heavy quark in doubly 
heavy systems. The symmetries and systematics we have discussed in the preceding sections should not be expected 
to work well when applied to baryons in which the strange quark is treated as heavy. Nevertheless, examining such a 
limit may provide clues to the understanding of the spectra of such baryons, particularly the Es. It may also provide 
some information on how badly the symmetries get broken as the masses of the heavy quarks decrease. 



1. Meson-Meson Ratios 



We begin this investigation with meson-meson ratios, eqn. (12). These ratios allow us to extract the effective 
energy scaling index, d, for the region of reduced masses involving strange diquarks. For strangeonia there is an 
appreciable amount of flavor mixing, and calculations of this mixing are inherently model dependent. Therefore, we 
rely on phenomenological models to estimate the mixing angles and to extract the masses of pure ss states. We treat 
three states as unmixed ss mesons: 0(1020), fta(1380), and the /2(1525). For the 77-77' mixing we use the mixing angle 
of 42° due to Thomas [57]; for the 1 3 P s $ state we use the mass of 1682 ± 4 MeV determined by Close and Zhao 
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|68j : and for the /i(1285)-/i(1420) mixing angle we use the value of 21° ± 5° due to a method by Close and Kirk 
[69j updated with more recent experimental data [70] . These mixing angles are similar to those found in a lattice 
calculation [70]. At this time, there are only enough experimental data to calculate the spin-averaged masses of the 
ground state and the first orbitally excited state of D s mesons. We calculate the orbital splittings of these mesons in 
table [X] The errors reported in these results only reflect the experimental uncertainties, and the theoretical error due 
to the choice of phenomenological descriptions is left undetermined. 



TABLE X: Spin averaged mass splittings for D s mesons and strangeonium from experiment [52]. The energy scaling rule, eqn. 
|tJ), predicts that the ratio of these splittings should be (jf L!L ) d We also extract the energy scaling parameter, d, from the 
ratio using the heavy quark masses from the RP model [3]. 



Splitting 


Am s¥ (MeV) 


Am c; (MeV) 


Am,;/Am,,; d 


M(1P) - M(IS) 


515.11 ± 7.23 


437.42 ± 0.56 


1.178 ± 0.017 0.617 ± 0.076 



The result for the scaling index d is close to the value that one would expect for a purely linear confining potential 
(d = |). This suggests that the nonperturbative nature of QCD has become more important for this region of reduced 
masses. We note that, unlike the ratios obtained in section VI A the value of d extracted depends on the quark model 
masses used, and the value of 0.617 is obtained using the quark masses of the RP model. This is significantly different 
from the scaling index for bottomonia and charmonia, where da 1 and the quark mass dependence vanishes. It is 
interesting to note that the orbital splitting of the cs system is nearly the same as the analogous splitting in cc and bb 
systems (see table [IT]). This suggests that dal scaling holds reasonably well for a range of reduced masses from that 
of the bb system down to that of the cs system, and then approaches d — | scaling appropriate to a linear potential 
for ss systems. 



2. Meson-Baryon Ratios 



With the value of d = 0.617 ± 0.076, we expect the ratio of the orbital splitting in the ss meson to that of the 
heavy diquark in a 5 baryon to be 1.534 ± 0.035. For the S c baryons, there are not enough experimental data to fully 
determine the orbital excitation energy of the sc diquark. The PDG [62] results show all three ground states, along 
with two low-lying negative parity states which might be interpreted as containing orbitally excited heavy diquarks. 
There should be five additional negative parity states associated with orbital heavy diquark excitations, but those 
states are not yet identifiable. For the S, there are only 2 low-lying, negative-parity states associated with an orbitally 
excited heavy diquark, and these two states are readily identifiable. Based on the RP model predictions, we expect the 
two lowest excited states to be orbital excitations of the heavy diquark; the S(1820) has quantum numbers J p = | , 

and we expect the S(1690) to have quantum numbers J p = \ based on general features of quark models. Treating 
these states as pure diquark excitations, we calculate the orbital splittings of the 5 baryon and the expected splitting 



derived from the meson-baryon ratios (eqn. 14) and the value of the scaling index extracted from strangeonium and 



D s meson orbital splittings (table DO. The predictions from this application of the supercolor symmetry are in table 



XI along with the spin-averaged prediction from the RP model. 



The experimental splitting and the splitting calculated using the supercolor symmetry are within 6% of each other. 
This agreement arises despite a number of possible issues with this treatment: (a) this prediction assumes non- 
relativistic kinematics for the quarks in the ss system; (b) the expected value for the splitting depends on the s and c 
quark masses used; (c) the solution to the energy scaling equation (eqn. [7]) from which we derived this result is only 
an approximation from which the value of d is chosen so that the solution reproduces experimental data for a finite 
range of reduced diquark masses; (d) we have ignored the mixing of heavy diquark excitations with excitations of the 
light degrees of freedom due to tensor interactions; and (e) the results for the 5 splitting depends on the models used 
to describe flavor mixing in strangeonium. It is also interesting to consider the prediction that one would obtain by 
assuming linear scaling (d = 2/3). In some sense this can be thought of as a non-relativistic, strong coupling limit. 
In this case the predicted splitting becomes Af-(lP) — M-(IS') = 324.5 ± 4.6 MeV, in slightly better agreement with 
the experimental result. 
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TABLE XI: S and H c heavy diquark orbital splittings. The Supercolor Symmetry predictions come from meson-baryon ratios 



(eqn. 14 1. We have used the estimate for d obtained from ss to cs meson- meson ratios (table [Xj) in rows 1 and 3 and the 
estimate for d from cs to cc meson-meson ratios in the last row. These estimates depend on the quark masses of the RP model. 
The second row uses a value of d = 2/3, which one would expect from a nonrelativistic treatment with a linear confining 
potential. 

Baryon Splitting d Supercolor Symmetry (MeV) RP Model Experiment (MeV) 

5 M(1P)-M(1S) °- 617± - 076 335.9 ±9.0 ^ 317 . ± 4 .7 

0.667 324.5 ± 4.6 

E c M(1P)-M(IS) 0^ ± 0.076 285.2 ± 6.5 ^ 

1.059 ± 0.004 209.9 ± 0.4 



3. Symmetry Multiplets 



The angular momentum multiplets obtained in the RP model are listed in tables |XH| and XIII for the flavor-sextet 
states and tables XIV and XV for the flavor-antitriplet states. The masses for the (H c , HQ and (3&, H'J states presented 
here differ from those published in ref. [3J, because the present spectrum has been recalculated using the SU(3)h e avy 
flavor basis, but the original work used the familiar SU (2t)u g ht basis. If (heavy) flavor-sextet and -antitriplet states 
are allowed to mix, the results would be identical with either choice of basis. The results for strange diquarks in 
the spinless model show a remarkable similarity to the multiplet structure of S cc , S;, c , and Hjj baryons. The energy 



splittings due to pairwise mixing (section V D I are slightly larger for baryons with strange diquarks than for diquarks 
comprised of b and c quark s. This type of mixing causes the broken degeneracies of the n^ Sd+ 1 (l ), n\ (l\) = 1 X P, 
IP multiplets in table |XII and the 1 3 P, IP multiplets in the "Spinless" column of tables XIV and |XV| The largest 



pairwise splittings are in the X l P\, 1-Pg multiplets. For strange systems this splitting is 41 MeV (for the S), and 
for doubly heavy systems the splitting is 24 MeV (for the S cc ). When spin interactions are included in the model, 
the largest spin splitting comes from the hyperfine splitting of the ground state multiplet. As the mass of the heavy 
diquark system is increased from the ss, cs, bs, cc, be, to the bb system, this splitting decreases from 195, 115, 87, 74, 
50, to 30 MeV, respectively. In the limit where the c and s quarks are infinitely massive the statistics of the heavy 
diquark should be irrelevant. Thus, the "E! c in the ground state antitriplet should be degenerate with the S c in the 
ground state sextet. In the spinless RP model, the mass difference between these two states is 1 MeV. In the full 
model, these states have masses that are within 3% of each other. 

Experimentally, only the ground state has been established. For the S c there are several excited states whose 
masses are known, but no quantum numbers have been assigned. It is therefore only possible to identify three of the 
angular momentum multiplets. The clearly identifiable E and S c multiplets are listed in the tables discussed in the 
previous paragraph, beneath the corresponding model prediction. For the 2 C ground state multiplet, we have identified 
the states which are primarily composed of flavor-sextet heavy diquarks by comparing masses with predictions from 
the RP model. Using baryon masses from experiment [55], the states belonging to the 5 ground state doublet are 
degenerate within 16%, while states belonging to the first orbitally excited multiplet are within 8% of each other. The 
S c ground state doublet is degenerate within 7%. 

At this point, it is possible to make a preliminary comparison between two different approaches for the S c (or S&) 
baryon: (1) treating the strange quark as a heavy quark and applying the heavy diquark symmetries that we have 
presented here; and (2) treating the strange quark as a light quark and applying HQS. As discussed above, treating 
the strange quark as a heavy quark leads to the prediction that the three lowest lying S c /S£. states will be degenerate. 
However, treating the strange quark as a light quark and applying HQS leads to the prediction that the three ground 
states will exist in a J p = (^ + ,| + ) doublet with J; — 1 (here, Ji refers to the spin of the light/strange degrees of 

freedom) and a J p — ( ^ + ) singlet with J; = 0. The properties of these two states are unrelated because each multiplet 
carries light degrees of freedom with different total angular momentum. To make this point more concrete, contrast 
the ground state of the S c /A c with the ground state of the E e /E' c . The spin averaged mass of the S c , J p = {\ + , | + ) 

HQS doublet is about 2500 MeV, the mass of the A c , J p = (\ + ) HQS singlet is about 2290 MeV, and the difference 
is 210 MeV. This difference in mass between the doublet (with J\ = 1) and the singlet (with J; = 0) illustrates the 
idea that these two multiplets are not related by any underlying symmetry. The analogous calculation for the S c /S^. 
baryons yields 150 MeV for the HQS doublet-singlet splitting suggesting the emergence of an additional symmetry. 
This additional degeneracy which HQS does not predict is easily understood in terms of the decoupling of the total 
angular momentum of the heavy-strange diquark. 
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TABLE XII: Symmetry multiplets for strange baryons in the RP model for flavor-sextet states with spin interactions removed. 
The decoupling of the diquark total angular momentum predicts that the states within a multiplet will be degenerate. 



Multiplet J p 




Mass (GeV) 




n2 d d+1 ( l p)jd> n A^A)j, (| Jd - Jt\, Jd + Jl) 









l 3 £i, 


ISi 

2 


/ 1 

(a 


4 ) 


(1.451, 1.451) 


(2.608, 2.608) 


(5.959, 5.959) 


1 1 D 


ISi 

2 


/ 1 

(a 


>2 ) 


(1.741, 1.741) 


(2.865, 2.865) 


(6.210, 6.210) 


2 bi, 


lb i 

2 


/ 1 
(a 


+ 3+\ 
'2 ) 


(1.9do, l.aoo) 


(3.064, 3.064) 


(6.404, 6.404) 


l s Di, 


ISi 

2 


/ 1 
(a 


+ 3+\ 
»3 ) 


(2.002, 2.002 ) 


(3.094, 3.094 ) 


(6.433, 6.433 ) 


l 3 D 2 , 


ISi 

2 


(1 


+ 5 + \ 
' 2 J 


(2.002 a , 2.002) 


(3.094 a , 3.095) 


(6.433°, 6.435) 


1 3 Do 


ID 1 

2 


(.a 


+ 7 + n, 
' 2 J 


1 Z.UUZ, Z.UUZ 1 


(3.095, 3.094) 


(6.435, 6.434) 


l 3 Si, 


lPl 

2 


(1 


- 3-S 
' 2 J 


(1.806, 1.806 6 ) 


(2.953, 2.953 6 ) 


(6.294, 6.294 6 ) 


l a S u 


IPs 

2 




3 - 5 - \ 
2 ' 2 J 


(1.806, 1.806 6 , 1.806) 


(2.953, 2.953 6 , 2.953) 


(6.294, 6.294\ 6.294) 


l 1 A, 


lPi 

2 


<i 

V 2 


+ 3+\ 

' 2 I 


(2.094 c , 2.078) 


(3.207 c , 3.194) 


(6.543 c , 6.531) 




IPs 

2 


( i + 

V2 ' 


3+ 5+\ 
2 ' 2 J 


(2.053 c , 2.094, 2.077) 


(3.174 c , 3.207, 3.192) 


(6.512 c , 6.543, 6.529) 


l 3 Si, 


2Si 

2 


a 
y 2 


+ 3+N 
' 2 J 


(2.091, 2.091) 


(3.226, 3.226) 


(6.560, 6.560) 


l 3 Si, 


ID 3 

2 


( i+ 

V2 ' 


3+ 5+\ 

a ' a > 


(2.132, 2.132, 2.129 d ) 


(3.268, 3.268, 3.262 d ) 


(6.602, 6.602, 6.596 d ) 


l 3 Si, 


IDs 

2 


( 3 + 
V2 ' 


5+ 7 + N 

2 ' a J 


(2.132, 2.129 d , 2.132) 


(3.268, 3.262 d , 3.268) (6.602, 6.596 d , 6.602) 



obcd These states mix significantly. 



4- Diquark Excitation and Flavor Symmetries 

The excitation symmetry emerges fairly clearly for diquarks comprising strange quarks. Table |XVI| shows the 
energy required to excite the light degrees of freedom in flavor-sextet DHBs in the spinless RP model f or mul tiplets 

shows 



XVII 



with two different diquark excitations: the ground state and the n/ d ^~ {l p )j d =l l Pi multiplets. Table 
the analogous splittings for flavor-antitriplet states. The excitation symmetry requires that the entries on the first 
row of table |XVI| should be the same as those on the second row, while those on the third row should be the same 
as those on the fourth row. In table |XVII| this symmetry means that entries on the first row should be repeated on 
the second, third and fourth rows, and entries on the fifth row should be repeated on the sixth, seventh and eighth 
rows. This symmetry should be broken by this particular model because of pair-wise interactions, but all of the 
excitation energies are within 17 MeV of each other for the flavor-sextet states, and within 25 MeV of each other for 
the antitriplet states. The effects of pair-wise forces are therefore seen to be small, even for diquarks involving the 
strange quark. It is also interesting to note that these relations work better for states which have identical quarks 
in the heavy diquark. This may be an indication that factorization emerges more readily for states containing such 
diquarks, and is a better approximation in those cases. 

Tables |XVI| and |XVII| also indicate the applicability of the heavy diquark flavor symmetry to strange diquarks. 
The excitation energies should be equal for any diquark flavor. Comparing the excitation energies of S baryons (table 



XVI) with those of baryons (table VII) reveals that the excitation energies are still remarkably similar. The worst 
disagreements occur for the Sf, baryons. Nevertheless, it is remarkable that these excitation energies are all so similar, 
as one would expect the 1/Md corrections to these symmetry predictions to be appreciable for the ss diquark. One 



other trend that can be seen in table XVI is that the diquark flavor symmetry works better when the masses of the 
quarks comprising the diquark are closer to each other. This trend is also apparent if we compare the splittings of 
E' bc baryons in table VIII with the splittings of E' c and E' b baryons in table XVII 
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TABLE XIII: Multiplets for strange baryons in the RP model for flavor-sextet states including spin dependent interactions. 
When available, the experimental masses [62] of these states are listed directly below the model predictions. The decoupling 
of the diquark total angular momentum predicts that the states within a multiplet will be approximately degenerate. 



Multiplet 


J p 




Mass (GeV) 




1 (l P )j d , n\(h)jt 


(| Jd — Jl\, Jd + Jl) 




1 — 'C 




l 3 Si, ISi 


V 2 '2 J 


(1.325, 1.520) 


(2.537, 2.652) 


(5.903, 5.990) 


2 




(1.318, 1.533) 


(2.469, 2.646) 




l 1 p l is I 

' 2 


/l- 3-N 
^2 '2 J 


CI 725 1 75Q1 


(9 S61 2 S60~l 


(H 207 6 206'! 






(1.690, 1.823) 






2 3 5i, 15 1 

2 


V 2 '2 J 


(1.891, 2.021) 


(3.030, 3.096) 


(6.377, 6.424) 


l 3 P>i, 15i 

2 


,1 + 3+N 
^2 ' 2 ^ 


(2.055, 1.934°) 


(3.119, 3.053") 


(6.450, 6.408°) 


l 3 D 2 , 15i 

2 


/•3+ 5+\ 
^2 '2 J 


(2.032°, 2.025) 


(3.113 a , 3.109) 


(6.446°, 6.443) 


1 3 D 3 , 15i 

2 


(5+ 2 + ^, 


(1.936, 2.035) 


(3.053, 3.115) 


(6.409, 6.447) 


l 3 5i, lPi 


\2 ' 2 > 


(1.826 c , 1.826 d ) 


(2.976, 2.952 d ) 


(6.311, 6.287°') 


l 3 5i, 1P| 


V2 '2 '2 ' 


(1.811°, 1.849 d , 1.883) 


(2.952, 2.963 d , 2.970) 


(6.287, 6.302 d , 6.307) 


^Pi, lPi 


n+ 3+n 

^2 ' 2 ^ 


(2.108°, 2.106) 


(3.206 6 , 3.195) 


(6.523 6 , 6.531) 


^Pi, 1P 3 

2 


(1 + 3+ 5+\ 
V 2 '2 '2 / 


(2.129° e , 2.115, 2. 137^ (3.162 DC , 3.206, 3.195) 


(6.508 6 , 6.542, 6.530) 


l 3 5i, 25i 

2 


\2 ' 2 ) 


(2.014 c , 2.144) 


(3.202 e , 3.256) 


(6.538, 6.579) 


l 3 5i, 1L» 3 

2 


(1+ 3+ 5+\ 
1 2 '2 '2 / 


(2.185, 2.149, 2.143 7 ) 


(3.277, 3.274, 3.267) 


(6.609, 6.606, 6.596) 


l 3 5l, LD 5 

2 


(-3+ 5+ 7+\ 
V 2 '2 '2 / 


(2.142, 2. 1480 2.148) 


(3.268, 3.270, 3.275) 


(6.598, 6.600, 6.607) 



These states mix significantly. 



5. Superflavor Symmetry 



In section VI F we used the superflavor symmetry to estimate that the energy required to excite the light degrees of 



freedom of a DHB from the ground state to the JJ 1 ' = | excited state is about 400 MeV. However, this estimate is 
only valid for very heavy diquarks. For strange diquarks, one can estimate the diquark mass as Md ~ Mq 1 + Mq 2 . In 
the RP model this gives a mass of about 1.1 GeV for the ss diquark. As such, corrections due to its finite mass could be 
quite important. This is illustrated in fig. [2] which shows the energy of orbital excitations for singly heavy mesons from 
experiment plotted against the inverse heavy quark masses taken from the RP model. These three data points come 
remarkably close to lying on a straight line considering the relatively small mass of the strange quark. These orbital 
splittings may be parametrized up to second order in the heavy quark expansion as AM = K 2 /Mq + K\/Mq + AM , 
where K\ and K 2 are constants which are chosen to reflect experimental results and AMo is the value of the orbital 
splitting that emerges in the infinitely massive heavy quark limit. Fitting the heavy meson data from table [IX| gives, 
K 2 = -4.26 x 10 7 MeV 3 , K Y = 1.80 x 10 5 MeV 2 , and AM = 388.9 MeV. HQET arguments indicate that K Y should 
be on the order of Aq C£) ~ 10 5 MeV 2 and K 2 should be on the order of A^ CD ~ 10 7 MeV 3 . Thus the values of the 
correction coefficients which we have obtained in this phenomcnological approach seem to indicate that inclusion of 
corrections up to 0{\/Mq) would describe these systems quite well. Because these corrections cannot be neglected 
in singly heavy meson systems one may speculate that they are also important for systems with strange diquarks. A 
systematic treatment of these corrections is left for possible future work. 



VII. CONCLUSIONS 



We have examined the systematics and symmetries of the spectra of baryons containing two heavy quarks, partic- 
ularly as they pertain to the model of Roberts and Pervin [3]- We have argued for an extension of the superflavor 
symmetry, which was proposed by Savage and Wise pQ, to include excitations of the heavy diquark, and we have 
examined some of the consequences of this extension for the spectra of DHBs. In this extension, the total angular 
momentum of the heavy diquark decouples from the light degrees of freedom independent of the excitation state of 
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TABLE XIV: Symmetry multiplets for the RP model for ~S! C states. When available, the experimental masses [62] of these 
states are listed directly below the model predictions. The decoupling of the diquark total angular momentum predicts that 
the states within a multiplet will be approximately degenerate. 



Multiplet 




1 


Mass (GeV) 


(lp)j d , n\{h)Ji 


/ 1 7 T| 7 
(\Jd — Jl\, Jd T 


- Ji) Spinless 


Spin 


1 D[) i ID 1 


(-2 J 






2 






(2.577) 


1 3 P , IS i 

2 


(D 


(2.865°) 


(2.907°) 


l 3 Pi, lSi 

2 


(r.n 


(2.865 a , 2.865) 


(2.814 a , 2.899) 


1 3 P 2 , lSi 


(f-.n 


(2.865, 2.863) 


(2.816, 2.890) 


2 1 5 , IS i 


(D 


(3.065) 


(3.032) 


l 1 ^, ISi 


/3+ 5+\ 
^2 ' 2 1 


(3.095, 3.095) 


(3.091, 3.090) 


l l So, lPx 

2 


(D 


(2.953) 


(2.927) 


1^0, IPs 

2 


(I") 


(2.953) 


(2.927) 


l 3 Po, lPi 

2 




(3.192°) 


(3.204°) 


l 3 Po, lPa 

2 


(l + ) 


(3.205 c ) 


(3.202 c ) 


I 3 Pi, lPi 


^2 ' 2 ^ 


(3.205 6 , 3.192 c ) 


(3.216 6 , 3.195 c ) 


I 3 Pi, 1P| 


(1 + 3+ 5+\ 
^2 '2 '2 J 


(3.205 6 , 3.192 c , 3.192 d 


) (3.212, 3.178 c , 3.205 d ) 


1 3 P 2 , lPi 


/•3+ 5+\ 
\2 1 2 ) 


(3.173 c , 3.192 d ) 


(3.209, 3.177 d ) 


1 3 P 2 , IPs 

2 


/ 1 + 3+ 5+ 7" 
V2 '2 '2 '2 


") (3.173 6 , 3.205 c , 3.205 d , 3.192) (3.116°, 3.206 c , 3.199 d , 3.S 


l^o, 25 1 

2 


(| + ) 


(3.221) 


(3.219) 


l^o, IDs 

2 


(| + ) 


(3.262) 


(3.241) 


l^o, IDs 

2 


(| + ) 


(3.262) 


(3.239) 



abcd These states mix significantly. 

the heavy diquark. Using the Feynman-Hellman Theorem, we have derived an equation which describes how the 
excitation energy of a doubly heavy meson or a heavy diquark scales with the masses of the heavy quarks and the 
strength of the interaction. The spectra of bottomonium and charmonium indicate that the excitation energy of these 
heavy-heavy systems are independent of the heavy quark mass, and grow linearly with the interaction strength. Wc 
have proposed a relation between excitations of heavy diquarks and those of doubly heavy mesons. This relation 
depends on the form of the confining forces. When combined with factorization, this leads to a relation between the 
spectra of DHBs and the spectra of doubly heavy mesons. Even though the model of ref. |3j was not constructed 
using any of the symmetries we discussed, we found that these symmetries nevertheless emerged. 

We have also explored the applicability of these symmetries to systems with strange quarks. Results from the RP 
model for strange diquarks are remarkably similar to results for DHBs. There is some indication that the consequences 
of factorization do not emerge as well for S c and baryons as they do for DHBs and the S. The fact that any of 
the consequences of factorization emerge is rather surprising and suggests that (at least for this particular model) a 
point-like diquark is not necessary for the system to factorize. S c and Hb baryons are also interesting, in that two 
different descriptions may be possible. On one hand, one may interpret their spectra in terms of the superflavor and 
supercolor symmetries that we have presented by considering these systems to be composed of a 'heavy' diquark. On 
the other hand, one may think of the strange quark as a light quark, in which case, the spectra can be interpreted in 
terms of the HQS and the SU(3)n g ht flavor symmetry. Whether these very different descriptions are complimentary 
or contradictory remains to be determined. Our estimate of the lowest orbital splitting of the 5 using the supercolor 
symmetry agreed with the experimental data. If this is not coincidence due to our choices of phenomenological models 
for strangeonium mixing, then it provides some evidence for the existence of this symmetry. Additionally, the orbital 
splittings of singly heavy mesons seem to show that HQS is applicable (phenomenologically at least) to strange quark 
systems. By extension, we expect that the superflavor symmetry will provide reasonable predictions for systems 
containing strange diquarks as well as heavy diquarks, but this remains to be tested. 

This work has concocted a description of DHBs in which the entire spectrum may be approximately constructed 
from the spectra of singly and doubly heavy mesons once the mass of the ground state DHB has been measured. We 



20 



TABLE XV: Symmetry multiplets for the RP model for S' b states. The decoupling of the diquark total angular momentum 
predicts that the states within a multiplet will be approximately degenerate. 



Multiplet 


J p 


1 — 'b 


Mass (GeV) 


n d Sd+1 ( l p)jd> nx(l\)j l 


(\Jd — Jl\, Jd 4 


- Ji) Spinless 


Spin 


^So, ISi 




(5.960) 


(5.944) 


l 3 Po, ISi 


(D 


(6.210°) 


(6.236°) 


l 3 Pi, ISi 


(I" 3~) 

V2 '2 / 


(6.210 a , 6.210) 


(6.173°, 6.228) 


1 3 P 2 , ISi 


(!".§") 


(6.210, 6.209) 


(6.173, 6.227) 


2 1 S , ISi 




(6.406) 


(6.390) 


UD2, ISi 


/ 3 + 5 + \ 
V 2 '2 / 


(6.435, 6.435) 


(6.434, 6.432) 


l l So, IPi 


(D 


(6.294) 


(6.285) 


I 1 So, 1P| 


(I") 


(6.294) 


(6.285) 


l 3 Po, IPi 




(6.529°) 


(6.540°) 


1 3 P , IPs 


(l + ) 


(6.541°) 


(6.539°) 


l 3 Pi, IPi 

2 


V 2 '2 J 


(6.542°, 6.529°) 


(6.547 6 , 6.531°) 


l 3 Pi, IPs 


(1 + 3+ 5+\ 
^2 '2 '2 J 


(6.541°, 6.529°, 6.529 d 


) (6.545, 6.514°, 6.536 d ) 


1 3 P 2 , IPi 


/•3+ 5+\ 
\2 ' 2 J 


(6.510 c , 6.529 11 ) 


(6.544, 6.513 d ) 


1 3 P 2 , IPs 

2 


/ 1 + 3+ 5+ 7" 
\2 '2 '2 '2 


") (6.510°, 6.540°, 6.541 d , 6.529) (6.462 6 , 6.542°, 6.541 d , 6.539) 


l^o, 2Si 

2 


(| + ) 


(6.555) 


(6.553) 


l 1 ^, 1-D| 


(| + ) 


(6.596) 


(6.587) 


USo, LD| 


(| + ) 


(6.596) 


(6.586) 



abcd These states mix significantly. 

have neglected spin-dependent interactions and assumed that factorization is valid. We have also ignored the mixing 
of heavy diquark flavor-singlets and diquark flavor-triplets due to hyperfine interactions. For the ratios involving 
doubly heavy mesons, we have also assumed that it is reasonable to use a non-relativistic potential model, or rather 
that such a model exists which captures the essence of QCD interactions for these systems. It would be interesting if 
ratios similar to the ones we have derived could be extracted from NRQCD. We have also assumed that the strength 
of interactions for a doubly heavy meson is exactly twice that of a heavy diquark. However, to our knowledge, there is 
no rigorous demonstration that this holds in general. It would also be interesting to see if there is a firmer conceptual 
foundation for the symmetry between doubly heavy mesons and heavy diquarks. Simplistically, this is a symmetry 
between a heavy quark and a heavy antiquark in a color singlet and two heavy quarks in a color antitriplet, but we 
have ignored how this relates to charge conjugation and parity inversion. The simple result that for bb, cc, and cs 
systems the energy scaling index d 1 reveals an important fact regarding doubly heavy hadron spectroscopy that 
models for baryons should not ignore - the energy required to orbitally excite two heavy quarks is nearly independent 
of the heavy quark masses, at least, for the range of masses that encompass these diquarks. Moreover, the heavy 
quarks available to hadron physicists are not massive enough for the spectrum to behave in a Coulombic manner. We 
hope that this work has demonstrated the insights into QCD which the study of doubly heavy hadrons offers. We 
propose to examine other aspects of the phenomenology of heavy hadrons in the framework of the symmetries and 
systematics that we have discussed. 
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TABLE XVI: Energy required to excite the light degrees of freedom of a strange SU (3) heavy flavor-sextet baryon from the 
ground state, J™ 1 = | + , to the excited state listed. We have averaged over the masses in each multiplet. This excitation energy 
should be independent of the spin, parity, and flavor of the heavy diquark, so that the energies in the first row should be the 
same as the corresponding ones in the second row, and the entries in the third row should equal those in the fourth row. These 
excitation energies should also be independent of the flavor of the quarks that comprise the strange diquark, so that all entries 
on a particular row should be equal. This is seen not only for the first four rows, but also for the last three rows. 



Diquark state 




Upper 


Multiplet 


Lower Multiplet 


m* - 


- mo 


(MeV) 


n 2 d Sd+1 {l P )j d 


n\(l\)j l 


( 


m*) 


(m ) 










lPi 

2 


(1 


".§") 


(f,§ + ) 


355 


345 


335 


1*P1 


lPi 

2 


(1 


+ 3+\ 
' 2 I 


(1- 3~ ) 
V2 ' 2 / 


345 


336 


327 


l 3 ^ 


IP 3 

2 


(I" 5 

\2 ' 2 


' 2 > 


fl + 3+N 
V 2 '2 1 


355 


345 


335 


1 X P1 


IP 3 

2 


(1+ 3 
V2 ' 2 


+ 5+\ 
' 2 I 


(I- 3~ ) 
V2 ' 2 / 


334 


326 


318 



l 3 ^ 


ID 3 

2 


(1+ 3 + 
V2 ' 2 


6+) 
' 2 > 


,1 + 3+N 
V 2 '2 J 


680 


658 


641 


l i S 1 


IDs 


(-3+ 5 + 
V2 ' 2 


' 2 > 


fl+ 3+N 
V 2 '2 V 


680 


658 


641 


l 3 ^ 


2Si 


(| + 


■1+) 
' 2 I 


fl+ 3+N 


640 


658 


601 



TABLE XVII: Energy required to excite the light degrees of freedom of a strange SU (3) heavy flavor-antitriplet baryon from 
the ground state, Jp = | + , to the excited state listed. We have averaged over the masses in each multiplet. This excitation 
energy should be independent of the spin, parity, and flavor of the heavy diquark, so that the energies in the rows 1, 2, 3, and 
4 should be equal, and the entries in rows 5, 6, 7, and 8. These excitation energies should also be independent of the flavor of 
the quarks that comprise the strange diquark, so that all entries on a particular row should be equal. This is seen not only for 
the first eight rows, but also for the last three rows. 
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